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PRELIMINARY INVESTIGATION OF AN ELECTRICAL 
NETWORK MODEL FOR ULTRASONIC SCATTERING 


James E. Maisel 
Cleveland State University 


ABSTRACT 

This study relates the behavior of acoustic attenuation in 
a solid to the electrical transmission line model where the elec- 
trical shunt conductance, which is frequency dependent, represents 
the loss due to the scattering sites in the solid. Results 
indicate that the absolute value of attenuation at a given fre- 
quency depends on both the normalized mean square deviation of 
the density and bulk modulus of the scattering sites from the 
ambient medium and the spatial scattering correlation function. 
Besides establishing the absolute value of attenuation, the 
spatial correlation function determines the attenuation profile 
as a function of frequency. 
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INTRODUCTION 


This report contains a study to simulate acoustic 
attenuation in a solid due to scattering by employing a one- 
dimensional electrical transmission line model. The electrical 
transmission model offers the advantage of studying attenuation 
as a function of frequency and scattering site density without 
the tedious process of preparing specimens. 

By judiciously assuming that the scattering mechanism is 
isotropic and the scattering sites are randomly spatially dis- 
tributed, it will be shown that the attenuation in a solid can 
be represented by a distributed frequency sensitive shunt con- 
ductance in the transmission line model. The density and bulk 
modulus of the solid will be represented by a distributed series 
inductance and shunt capacitance, respectively. 

Since this study is attempting to relate a three-dimensional 
phenomena, scattering in a solid, to a one -dimensional model, 
electrical transmission line model, the concept of effective 
scattering cross section will be very important in the analysis. 
The collapsing of a multi -dimensional problem into a lower 
dimension simplifies the complexity of analysis providing that 
the desired salient features are retained. This is exactly 
what is done when measuring attenuation due to scattering in 
a physical specimen using a one -dimensional acoustic wave. 
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Assuming the above assumptions are valid, the concept of 
normalized mean square deviations of density and bulk modulus 
of the scatters from the ambient or host medium and the spatial 
correlation function of the scattering sites will be shown to 
be very important in describing the attenuation as a function 
of frequency. Results will predict Rayleigh and viscous type 
absorption. 

The absolute level of attenuation will depend on the 
normalized mean square deviation of the density and bulk modulus 
as well as the spatial correlation function; while the profile 
of the attenuation as a function of frequency will depend 
strongly on the Fourier transform of the correlation function. 
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I . ACOUSTIC WAVES 


Propagation of an acoustic wave is governed by two 
fundamental equations : the equation of motion and the equation 
of conservation of mass. The equation of motion is 

p a? +Vp = 0 a-D 

where p is mass density, v is velocity and p is pressure. Because 
the total time derivative is used in Equation (1-1) , it is the 
time rate of change in a coordinate system moving with a parti- 
cular portion of the medium. This is called the Lagrangian 
description. The Eulerian description is the partial time 

g 

derivative ^ at a "fixed point" in space as the medium moves 
past that point. The total and partial derivatives are related 
by the following operator expression 

= Tt + (1_2) 

The variables p, p, and v can be expressed in terms of 
their average values, p Q , p Q , and v Q and the small vibrating 
acoustic wave components p^ , p^, and v^. Assuming the medium 
is stationary (v o =0) and the magnitudes of the small vibrating 
acoustic wave components, p-^ and p-^ , are small compared with the 
average values p Q and p Q , Equation (1-1) becomes 

3v, 

p oTT +v Pl = 0 d' 3 ) 
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where S-Vv = (v Q +V * 1 * * * V (v 0 ^) * V V <v 0 > * 0 because the 

medium is stationary. 

The conservation of mass is given by 

V - (p$) + || = 0 (1-4) 

where V • (pv) = v- Vp + pV • v= p Q V ■ v^ and 3p/St = Sp-^/St. 

Hence, Equation (1-4) can be rewritten 

p o V-v x + 3p x /3t = 0. (1-5) 

In general, pressure p is a function of density P. For 
small p-^, it can be expressed by the following relationship 

P = P Q + Pi = p o + (3P/3p) p Pi- ( 1_6 > 

2 2 

Thus p-^ = (3p/3p)p p 1 = c P 1 where c = (3p/3p) p . The value 
of c is the velocity of the acoustic wave in the medium. Note 
if the density is independent of pressure (3p/3p=0) the velocity 
of propagation is infinite. Equations (1-3) and (1-5) can 
be written in terms of two unknowns: p-^ and v-^ . 

3v, 

p o TE~ + 7 Pi = 0 (1-7-a) 


1 3P 1 

^--3T + V - v l = ° 

c p o 


(1-7-b) 


Taking the divergence of v^ in Equation (1-7-a) and substituting 

in (1-7-b) it can be shown that 


i LZl - v(^i) =o. 


c Po at' 


Po 


( 1 - 8 ) 
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Equation (1-8) is the basic acoustic wave equation for a 
stationary medium. If the medium is uniform (p Q = constant). 
Equation (1-8) can be simplified as 


1 3 2 Pl 
-J 


- V' 


V 1 = 0. 


(1-9) 


c“ St*" 

For a time harmonic case [exp(-jcot)] , Equation (1-9) becomes 


2 

a_> 

X 


(v z + ^)P X = o 


d-io) 


where P^ represents the phasor pressure and is a function of 
space and co the angular frequency (co =2irf) . The phasor velocity 
can be determined from Equation (1-7-a) 

- VP 

V, = T-i-. d-11) 

1 Jo)p 0 

-V 

The power flux density vector S is given by 


S = 2-P x V* (1-12) 

where * represents the conjugate operator. 

For a medium that is infinite in extent, a pressure wave 
traveling in the +x direction is given by the following ex- 
pression 

P 1 = P 1q exp [ Yx ] (1-13) 

where P-^ o represents the amplitude of the pressure wave and 
Y = a+ jk, the propagation constant associated with the pres- 
sure wave. The real part of Y (a = Real [Y] ) measures the 
attenuation per unit length while imaginary part of y (k=Im [y] ) 
measures the phase shift per unit length. 
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The amplitude of the velocity, V^ Q can be expressed in 
terms of P^ Q by substituting Equation (1-13) into Equation 
(1-11) . The result is 



* YP lo 

3<up 0 


(1-14) 


where x represents a unit vector in the x-direction. 


The characteristic impedance Z Q of the medium is defined 
as the ratio of P-^ o /V^ o and in general will be a complex ex- 
pression. The characteristic impedance is 



(1-15) 


If the medium is lossless (a=0) , the characteristic acoustic 
impedance will be a real expression (Z Q = jeup 0 /jk = P Q c) . 
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II . ACOUSTIC SCATTERING 


Consider a particle, as shown in Figure l-a, illuminated 
by an incident acoustic plane wave propagating in the x-direction 

P i (x) = P io exp(jkx) (2-1) 

where k = w/c. Assume in the discussion that follows that 

P io - X - 

The scattered acoustic pressure wave is given by 

P .<?•) = f ( 1 , 0 ) ex P (i kr > for R>^~ (2-2) 

where A is the wavelength in the medium; D is the particle 
dimension; f(i,o) represents the information about the ampli- 
tude, phase and polarization of the scattered wave in the far 

2 a /v 

field (R > D /A) ; and i and o are unit vectors associated with 
the incident and scattered wave respectively. [1 ]. 

The incident and the scattered power flux in a lossless 
medium are given by 

S. = i(|P.l 2 / 2 Z o) (2-3-a) 

S s = S(|p s | 2 /2Z o ) (2-3-b) 

where Z q is assumed to be a real quantity. It will be shown 
later that for a low loss medium, Z Q is essentially a real 
quantity . 

The differential scattering cross section is defined as 
follows , 
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(2-4) 


a d (i,o) = lim[R 2 S /S.] *= lf(i,o)I 

u r’-^-oo 


The above equation represents an effective geometrical area 
that the incident wave "sees" and this area is equal to the 
magnitude squared of the scattering function f(i,o). A more 
meaningful definition of a cross section is one that accounts 
for the scattered power at all angles (in all directions) 
surrounding the scattering site. This cross section is de- 
fined as the scattering cross section a and is given by 

s 


= fa* dfi = f ' f(l,o) 


d n 


(2-5) 


4-rr u 4ir 

where d^ is the differential solid angle. 

If both scattering and absorption are present, the sum 
of the scattering and the adsorption cross sections is 


+ a. 


( 2 - 6 ) 


where a & represents the absorption cross section and is equal 
to the volume integral of the loss inside the scatterer: 


o = f WdV/S. 

a V 1 


(2-7) 


where W represents the power density loss inside the scattering 
site. 
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HI. THE FUNCTIONAL RELATIONSHIP OF ° s WITH RESPECT TO 
WAVELENGTH AND THE VOLUME OF THE SCATTERING SITE 


According to Equation (2-2) the amplitude of the far 

field pressure wave is inversely proportional to the distance 

R. This must be the case because the total power radiated from 

the scattering site must be independent of radius R (area of 
2 

a sphere = 4 ttR ) . 

The scattered pressure wave P g is a linear function of 

the incident pressure wave P. and thus P at a distance R is 
r is 

proportional to P^ and the volume V of the scatterer: 

[P s |>| P .| liitM = |p ± | constant V . (3-1) 

-2 

The constant in Equation (3-1) should be of dimension L , and 

since it is a function of wavelength, it should be proportional 

_ o 

to (wavelength) . Hence, 

o s ~ jf(i,o)| 2 ~ V?x 4 . (3-2) 

The above analysis assumes that the size of the scatterer is 
smaller than a wavelength and is referred to as Rayleigh scat- 
tering . 

In order to determine if the geometry of the scattering 
site has an effect on the results as stated in Equation (3-2) , 
let the scattering site be a cylinder whose diameter is small 
compared to X and whose length approaches infinity. (Obviously 
there is no such site in the "physical world.") The amplitude 
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II 


of the far field pressure wave must be inversely proportional 
to vH because the total power radiated per unit length of the 
cylinder must be independent of radius R (area of cylinder per 
unit length = 2irR) . Hence 

A A 

IP 1 = |P. ! l « |p.| (constant) (^2.^ ( 3 - 3 ) 

3 1 1 Jr 


where a is the radius of the cylinder. The constant in 

-3 IT 

Equation (3-3) should be of dimension L ' , and since it is 

a function of wavelength, it should be proportional to 
-3/2 

(wavelength) ' . Hence, 


a g ~ | f (i,o) 


(cross sectional area of cylinder)' 

73 


(3-4) 


For very long wavelengths and first-order approximation, 
Morse and Ingard [2] have shown that the scattered power flux 
is 


K- 


f 3 a 4 


S i (l-2 cos 0 )' 


(3-5) 


where K is a constant of proportionality and 0 is the angular 
displacement . 

The results of the above analysis indicate that the 
scattering cross section is proportional to the frequency to 
the fourth power for a site whose dimensions are small compared 
to wavelength and the third power for a site that has one 
dimension (e.g. the length of the cylinder) approaching infinity. 
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IV. SCATTERING FROM AN INHOMOGENEITY 


The energy lost by the incident wave may be absorbed or 
scattered by a scattering site. In Section II, it was shown 
(Equations (2-5) and (2-6)) that the amount of energy lost 
per second by an incident plane wave divided by the incident 
wave's intensity or power flux is called the cross section a. 
of the object. If there are, on the average, N scatterers 
per unit volume of the medium, the incident wave, for each unit 
volume of travel, looses Na t of the wave's power flux. Hence 
the intensity of a plane wave propagating through a scatterer- 
filled medium is 

I Si = JSJ exp(-Nd t x) (4-1) 

where x is the distance of penetration. The attenuation 
factor is ^-N a t - 

Morse and Ingard [3] have investigated the propagation of 
a pressure wave in an inhomogeneous medium. The essence of 
their investigation will be presented. Restricting the 
analysis to a single frequency (co/2tt) and supposing there is 
a region "R" , which can be inscribed in a sphere of radius "a" , 
the density and compressibility for the region inside and 
outside the sphere can be described by a set of quantities 

(p e , I?! £ R 

Density = ^ | ? | > R (4-2-a) 
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- We.c.*) 


Compressibility = 


K = l/(pc 2 ) 


! r < R 


r| > R 


(4-2-b) 


where c and c and P and p are the wave velocities and mass den- 
e e 

sities inside and outside the sphere, respectively. 

The time- independent equation for the pressure amplitude 


is 


V 2 P + k 2 P = -k 2 Y K (?) P - V- [Y p (?)VP] 


(4-3-a) 


where 


k = (o/c 


(4-3-b) 


y K 


i?i<r 

( 0 . |r|>R 


Y 


P 


P e -P 


1 

/ 


!?! <r 
\r\ >R. 


See Figure 1-b. 

At large distances from the sphere (r > > a, the radius of 
the sphere enclosing region "R”) , the pressure wave would be 

p = p. [exn(ik x)+f (£ , o) exp (jkr) /r] C4-4') 

where the first term in the brackets represents a plane wave 
traveling in the x-direction while the second term represents 
a scattered wave. The key term in Equation (4-4) is f(i,o), 
the angle-distribution factor, and is given by the following 
expression , 
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f(i,o) = pr J (Y K P-j(Y p ) £ • V'P)exp(-jk6.r') dv' . (4-5) 




If a small particle is assumed, the first term inside the 
integral in Equation (4-5) gives isotropic scattering and the 
second term gives scattering proportional to cosine of the 

A A 

angle between the unit vectors i and o. 

Unfortunately, Equation (4-5) depends on the total pres- 
sure wave (incident and scattered pressure waves) and therefore 
Equation (4-5) is not a complete description of the scattering 
amplitude in terms of known quantities. 

Chernov [4] analyzes the problem by the method of small 
perturbations. Assuming small changes in density and the 
velocity of propagation, the scattered pressure wave can be 
expressed by the following equation, 


P 

s 



f [2k 2 

V 



j k 


] r ex P(j k < r+x ') dv 

(4-6) 


where r is the distance from the scattering element (x'y'z') 
to the observation point x,y,z, P^ o represents the amplitude 
of the incident pressure wave, and Ac/c Q and Ap/p Q are the 
normalized change in the wave velocity and density. Asstiming 
Ap/p o is independent of position within the scattering site, 
the second term in the bracket in Equation (4-6) is zero. Since 
the average velocity c Q = /B Ao /P Q where p Q and B Aq are the 
average density and bulk modulus, respectively, the normalized 
velocity change is 


Ac _ 1 .Ajo _ ^ B A . 
c o ^ p o B Ao 


(4-7) 


where AC=C-C e , Ap = p -p e> AB^ = - B A , and 
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Ap/p o and AB A /B A( ^ represent the normalized change in the density 

bulk modulus. See Appendix A. Treating Ap and AB A as random 

flucuations and introducing the normalized correlation 

functions N p (r) and (r) where N p (r) = ftp^Ap^/Ap I 2 and 

.A 

N b (r) = AB a] A B a2 / AB a , and assuming zero correlation between 

density and bulk modulus, the mean-square scattering pressure 

P P* can be determined using Equation (4-6) . The result is 
s s 


1 o 1 




32TrR^ V X V 2 


:^2T n P< ? > + ^2^ N B a ]exp[jK*r] dv x dv 2 


AB, 


B 


Ao 


(4-8) 


Where — * — (R is the distance from the origin of coordinates 
r R 

to the observation point. See Figure 2. 
r = r^ - r^ the vector separation between two points 
inside volume V. 

K = k (n - n^) 

I K I = 2k sin 0/2 where 0 is the scattering angle. 

2 

Ap = mean- square density. 

2 

AB a = mean- square bulk modulus. 

In order to obtain a better understanding of Equation (4-8) , 

asstime that N p (r) = Ng (r) . Using transformations similar to 

A 

the types used by Chernov [4 ] , the result is 

I P S 1 2 = A I P io I 2 \ [ ^/Po + ^A /B Ao ] f N(r) sin (kr) rdr (4-9) 
R 0 


where A is a constant and the upper limit on the integral has 
been set to equal infinity, since the dimensions of the 
medium are assumed to be large compared to the correlation 
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distance. Note the mean-square pressure in Equation (4-9) is 
a function of the scattering angle 0 . 

The scattering area a g , which is the ratio of the total 
power scattered to the incident power density, can be de- 
termined from Equation (4-9) assuming a statistically isotropic 
medium and integrating the scattering angle © from 0 to ^ 
radians . The result is 


a 


s 


= Ak Z t Ap Z / p ; 


AB i/ B io ] 


( [l-cos(2kr) N(r)]dr 
0 


(4-10) 


In order to gain insight about the behavior of a as a 

y. S 

function of k, let N(r) = e -r ^ a and N(r) = e - ^ " where "a" 
can be interpreted as the "space constant" of the medium. The 
results are 


N(r) = e 


-[r/a| 


a = A [ 4p 
s 


Vp2 + 4 ^/b| 0 ] 


4a 3 k 4 

l+4a Z k Z 


(4-11-a) 


N(r) = e- ,r/a|2 : a s = At A / p 2 Q + AbJ/B Z q ] ( ^a)k 2 (l-e~ (ka) ) 

(4-11-b) 

See Appendix B. 

For ka<<l the scattering area is proportional to k 4 and for 

o 

ka >:> l the scattering area is proportional to k for both 
correlation functions . Figure 3 is a plot of the following 
two functions which are the frequency dependent portion of 
Equation (4-11) . 


f,(ka) = (ka) 

1 l+(2ka) 2 


f 2 (ka) = (ka) 2 (l-e' (ka)2 ), 


(4-12-a) 


(4-12-b) 


16 



Comparing the slopes of f^(ka) and f ^ (ka) with the function 
(ka) n where n = 2, 3 and 4, it can be seen that there is a 
gradual change in n as ka varies ; Rayleigh for ka << 1 (n=4) 
stochastic scattering for ka>>l (n=2) and a transition region 
of ka ~1. 
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V. SIMULATING SCATTERING USING 


A ONE -DIMENSIONAL TRANSMISSION LINE 

Propagation of longitudinal sound waves in a solid and 
the propagation of transverse electromagnetic waves on an 
electric transmission line obey the same second order partial 
differential equation assuming the amplitude of the sound 
wave is small. The two parameters that characterize the medium 
in which the wave exists will be attenuation per unit length 
(a) and phase shift per unit length (k) . These two parameters 
can be combined together using complex notation to generate 
a term called the propagation constant of the medium which is 
expressed as 

Y = a + j k (5-1) 

where k = 2^/^ = 2irf /c. 

A low-loss medium is one that absorbs very little energy 
from the wave as it propagates through the solid and can be 
mathematically defined when a<< k. For metals, 0.1 <_10 tiepers/cm 

and k ~ 625 radians/cm at f=50MHz. Thus metals can be classified 
as essentially a low-loss material. 

Before the model of a wave traveling in a low loss medium 
is discussed, a short discussion on low-loss transmission lines 
will be presented. 

Figure 4 illustrates an electrical T-network that 
is terminated in its characteristic impedance, where the 
acoustic impedance pc will be characterized by analogy with 
the electrical impedance Z Q (Z o =pc) . It can be 
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shown [5 1 that when the circuit is terminated in its 
characteristic impedance, the input impedance is equal to 
the characteristic impedance which in turn is a function of 
the per unit impedance Z, per unit admittance Y and the 
physical length Ax. The characteristic impedance is 


Z o " 


# <1 + 


ZY 


Ax ) 


2 .. 1 % 


(5-2) 


The propagation constant can be expressed in terms of Z 
and Y as 

e YAx = e ctAx e jkAx = 1 + /gy Ax + j-(ZY)AX 2 + (5-3) 


From Equations (5-2) and (5-3), it is obvious that the 
characteristic impedance and the propagation constant depend 
on the per unit series impedance , per unit shunt admittance 
and the length of the transmission line section in a rather 

complicated manner. If the low-loss constraint is imposed by de- 
leating the series per unit resistance, then the per unit acoustic 
impedance and admittance are: 

Z -.jup o (5-4-a) 


Y = G + j*VB Ao (5-4-b) 

where P Q , B AcV w and G re P rese nt the density of the material, 
the bulk modulus, angular frequency (2irf) and the per unit 
conductance, respectively. If G^co/B^and yAx << 1 , then 

Z o"v/P 0 R Ao (5-5-a) 

and 

P /pi 

(r Ax)=exp (aAx) exp (jk A x) ~1 + >7 / p o B Ao Ax + jco /-2 ax. (5-5-b) 

^ d Ao 
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The left side of Equation (5-5-b) can be approximated by 
1 +TAx = 1 + aAx + jkAx assumingYAxccl . The attenuation and 
phase shift per unit length is 


and 



(5-6-a) 



(5-6-b) 


The results shown in Equation (5-6-a) indicate that a is 
linearly related to G and k is equal to <o/c where c = 
From Section IV, it was shown that 


a 


ht a . 
2 s 


(5-7) 


Equating Equations (5-6-a) and (5-7) and solving for G, the 
result is 


G = 


/p o B Ao 


(5-8) 


The conductance G can be expressed in terms of Equation (4-11) . 

NA 


N(r) = e“ r/a 


: G = 


JPn* 


o Ao 


Ap 


•D 2 

a b a 


2 + ZZ 


B 


Ao 


4a 3 k 4 


l+4a 2 k 2 


(5-9-a) 


N (r) 






(1-e 


- (ka) 


2 

). 


(5-9-b) 


The conductance in Equation (5-9) reflect the scattering loss 
only and not absorption loss. As the mean-square values of 
Ap 2 and AbJ approach zero (perfect medium) the conductance 
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and attenuation become zero (no losses due to scattering) . 
Equation (5-9) also indicates that the conductance will be a 
function of frequency and the exact relationship will depend 
on the correlation function. From Figure 3, it appears for the 
two correlation functions assumed, the behavior of the con- 
ductance as a function of frequency (ka = 2irfa/c) would be 
approximately the same. 
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VI . CALCULATION OF ATTENUATION VERSUS FREQUENCY 


USING THE ONE-DIMENSIONAL TRANSMISSION 
LINE MODEL 

A computer program was written based on the electrical cir- 
cuit model shown in Figure 4. Table I lists the parameters used 
in the program and their corresponding numerical values. These 
values represent a medium having the approximate characteristics 
of steel. 


TABLE I 



• 

Symbol 


Density 

P 

7.84' 10 3 kg/m 3 

Young ' s modulus 

b a 

19. 6* 10^0 newton/m 3 

Length of 
Transmission line 

Ax 

10 B meter 

Velocity of 
Propagation 

V 

3 

5 '10 meters/sec. 

Average separation 
between scatterers 

a l 

a 2 

5'10 - ^ meter 
15 •lO - ^ meter 

— — l 


Two average correlation spacings (a-^^) were chosen in 
order to investigate the effect that superimposing two media 
had upon each other. The numerical choice was based primarily 
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on having the attenuation exhibit the Rayleigh and stochastic 
scattering in the 20 to 100MHz frequency range. The strength of 
each correlation spacing was controlled by weighting factors 
and W 2 which are associated with a^ and a 2 » respectively. 

The weighting factors and W 2 can be introduced in the 
electrical transmission line model by replacing the conductance 
G in Figure 4 by two parallel conductances such that G = W^G^ + 
^2^2 w ^ ere Gq and ^2 account for ea ch set of scatterers associated 
with a^ and a 2 respectively. By controlling the values of 
and W 2 (with the constraint that W-^ + W 2 = 1 for all cases) , a 
blend of each set of scatterers can be obtained. 

In order to circumvent the statistical determination of the 
structure of the medium (for example the variance of density and 
Young’s modulus A and AB^ , respectively), it was assumed that a 
medium with a-^ = 5 microns has an attenuation of 1 neper/cm at 
a frequency of 20 MHz. This assumption is based on the results 
of Vary [6], The variances Ap^ and AB^ were assumed to be in- 
dependent of the correlation spacing. This was based on the 
fact that the variances act as a scaling factor for attenuation, 
while the correlation function contributes to the attenuation 
profile when the behavior of attenuation is studied as a function 
of frequency. 

Figure 5 represents attenuation versus frequency for the 
two average separations in Table I. The medium with the higher 
scattering density (a^ = 5 microns) exhibits the higher 
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attenuation at any frequency. The general attenuation profile 
for both correlation separation is consistent with theory. At 
low frequency the exponent of the frequency is slightly less 
than 4 (Rayleigh) and greater than 2 (stochastic scattering) at high 
frequencies. Tripling the average correlation distance reflects 
a decrease in the scatterer density and a change in the attenuation 
profile. This is consistent with the results shown in Figure 5. 

The effect on attenuation of combining two media via 
weighting factors and W 2 are shown in Figure 6 . Three dif- 
ferent frequencies were chosen to study the effect different 
weighting factors have on attenuation. All three curves exhibit 
a negative slope which is consistent with changing from a medium 
of high scattering density to a medium with a low scattering 
density. Defining the attenuation-weighting factor sensitivity 
as a change in attenuation per change in weighting factor, one 
immediately sees, according to Figure 6, that the largest nega- 
tive sensitivity occurs at f = 100 MHz and when the weighting 
factor is zero and W 2 i- s unity. 
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VII . DISCUSSION 


An electrical transmission line model of a one-dimensional 
acoustic wave propagating in a solid has been developed to 
study the behavior of attenuation as a function of frequency . 

The equivalence of energy loss due to scatterers in the solid 
can be modeled as a conductance that is dependent on frequency . 
This conductive element causes the electrical wave to experience 
attenuation as it propagates along the transmission line which 
is similar to the phenomenon that occurs when the acoustic wave 
propagates in a solid. 

Theoretical results developed in this report indicate that 
if the distribution of the scattering sites is random, attenuation 
as a function of frequency depends on the functional behavior 
of the spatial correlation function and the normalized mean 
square deviation of the density and bulk modulus of the scat- 
ters from the ambient solid. The absolute level of the 
attenuation is primarily determined by the sum of normalized 
mean square deviation of both density and bulk modulus of the 
scatterers and the spatial correlation function, while the 
ptofile of the attenuation versus frequency is determined 
principally by the spatial correlation function. Assuming 
correlation functions of the form e“^ r/a ' and e~\ r/a ^ 2 , the 
results indicate that attenuation depends on the fourth power 
of frequency if ka<< 1 where k is the wave number and "a" 
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is the average separation of the scattering sites; while for 
ka >> 1 the attenuation depends on the second power of fre- 
quency. Between these two extremes, attenuation is a rather 
complex function of frequency . 

It was shown that the attenuation, which is directly 
proportional to the net cross sectional area of the scattering 
site, is proportional to both the square of the wave number k 
and the difference between the Fourier transforms of the scat- 
tering site correlation function. The Fourier transforms of 
the correlation function that form this difference are evaluated 
at k = 0 and k f 0, respectively. Results indicate that for 

the Fourier transform evaluated at ka>>l, the attenuation is 

2 

due to viscous type absorption ( a a f ) . For the case when 
ka<<l, the absorption is Rayleigh (a^f^) . In essence, the type 
of absorption mechanism depends strongly .on frequency components 
that are contained in the correlation function. 

The results of this investigation were tested by assuming 
a medium with a base level of 1 neper/cm at a frequency of 20MHz 
and a = 5 microns. The spatial correlation function, e~^ r ^ a ^, 
was assumed in all cases. Based on the above attenuation con- 
straint, attenuation curves were generated using the electrical 
transmission line model for a = 5 and 15 microns over a fre- 
quency range from 20 to 100 MHz . The attenuation profiles 
developed follow the typical measured attenuation of a physical 
medium. For a = 15 micron medium, the behavior of attenuation 
as a function of frequency was less than for the 5 micron 
medium. This is consistent with a physical medium since "a" 
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represents the average separation between scattering sites . An 
increase in "a" represents a decrease in scattering site density. 

Two media, a = 5 and 15 microns, were superimposed with the 
mixture controlled by weighting factors and The weighting 

factors W t and are associated with media a = 5 and 15 microns, 
respectively. By varying and W 2 with the sum of and W 2 
equaled to unity, the blend of the media can be varied. At- 
tenuation versus weighting factors was investigated for fre- 
quencies of 30, 50 and 100 MHz. Results indicate that the 
greatest attenuation sensitivity occurs at f = 100 MHz , -+■ 0 

and W 2 1 . This is consistent with the physical medium, since 
the greatest change in attenuation occurs when the scattering 
site density is slightly greater than zero and the frequency 
is at its upper bound (f = 100MHz in this investigation) . 

This study indicates that future efforts should focus on 
comparing the attenuation as determined by using the pulse- 
echo- technique with the above approach. Assuming the medium 
has single type scattering sites and the scattering is isotropic, 
the normalized mean square deviation of the density and bulk 
modulus scattering sites can be determined by identifying the 
density and bulk modulus a priori and the corresponding 
spatial ensemble of these sites. The spatial ensemble can be 
determined by studying functional behavior of density and bulk 
modulus along radial lines eminating from a point on the sur- 
face of the specimen. This assumes that the surface configuration 
of the scattering sites will be representative of the distri- 
bution of these sites inside the medium. Details will have to 
be investigated on exactly how this information is to be obtained. 
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Knowing the functional behavior of the spatial ensemble 
as described in the above paragraph, the spatial correlation 
function can be determined using techniques found in statis- 
tical analysis. In this investigation, two correlation functions 

that are commonly found in physical systems, e - ^ r / a ^ and 
-lr/al 2 

e , were chosen. It seems plausible that correlation 

functions similar in behavior to the above functions should 
occur, with modifications, in order to account for specific 
spatial distribution of scattering sites. 
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VIII . CONCLUSION 


It has been demonstrated that using the electrical trans- 
mission line model for studying the scattering mechanism by 
ultrasonic waves in a solid is feasible if a frequency sensitive 
shunt conductance is introduced into the electrical model . The 
presence of the shunt conductance causes the electrical signal 
to be attenuated as it propagates , which is similar to the be- 
havior of a ultrasonic wave propagating in a solid. 

Mathematical results indicate that if a solid medium is 
viewed as realization of an ergodic spatial random process, the 
attenuation is related to the Fourier transform of the correlation 
function . 
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APPENDIX A - MODEL OF HOST MEDIUM 


Consider the model of a medium shown in Figure A-l-ra. 

The host medium can he characterized by a mass density p and a 
Young’s modulus B^, while the scattering material by p g and B^ e . 
The sectional view shown in the Figure depicts the scattering 
sites by islands surrounded by the host medium. Only one hand 
of scattering material is assumed which implies that all islands 
are identical. 

Assuming the scattering sites are randomly distributed 
spatially, the statistical characteristics of the distribution 
should be independnet of the orientation of the slice through 
the medium. Also, it will be assumed that the realization 
shown in Figure A-l-a is one realization from an infinite num- 
ber of possibilities of a random process and that this process 
is stationary and ergodic to at least the second moment. This 
implies that ensemble and spatial averages are the same and the 
ensemble and spatial correlation functions are equivalent. 

In other words, all statistical information up to and including 
the second moment is contained in any specimen composed of the 
same host medium and the same randomly distributed scattering 
sites . 

An ensemble of realizations as a function of space can be 
generated by moving along a line through the composite medium 
such as is shown in Figure A-l-b. Here the mass density is 
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plotted as a function of space. Young's modulus could be 
studied in a similar manner. Because the random process was 
ass um ed to be ergodic to the second moment, just one of the 
realizations would be sufficient to determine the average, 
variance and correlation function for the random process pro- 
vided that the dimensions of the medium is large compared to 
the spatial fluctuations in mass density or Young's modulus. 

The following mathematical development will center on 
the behavior of the mass density in the specimen. Parallel 
results can be obtained for Young ' s modulus . 

There is an average density, p Q , associated with the 

random process and it is shown in Figure A-2. The variance 
— 2 * 

Ap measures the amount of spread p has about p 

In order to completely describe the process to the second 
moment, the correlation function must be introduced. It measures 
how well the density at one point in the specimen compares to 
another point. Since the process is assumed to be ergodic, it 
is the separation distance that is important and not the ab- 
solute distance as measured with respect to some* origin. The 
correlation function can be expressed in the following manner 

Ap-, Ap 2 = Ap 2 N p (r) (A-l) 

where Ap z and N (r) represents the variance and the normalized 

r 

correlation function, respectively. These concepts were used 
in Section IV. 
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APPENDIX B - DEVELOPMENT OF EQUATION (4-11-a) 


This appendix presents the development of Equation (4-11-a) 
starting with Equation (4-9). For pxirposes of simplification, 
let B = [ Ap^/p^ + • Thus Equation (4-9) can be ex- 

pressed in the following manner: 


'0 irR 


AiB , ® N.(r) sin (2Kr sin (S-) ) 

“71 k j " //\\ rdr 


2k sin (0/2) 


(B-l) 


where 0 is the scattering angle, a represents the cross 

s 0 

sectional area for a given scattering angle and equals A/2 
in Equation (4-9) . The net cross sectional area is determined 
by integrating Equation (A-l) over a sphere of radius R. Hence 


7r 9 

a - Co 2ttR sin 0d0 = 
s o s 0 


iB o « it sin(2kr sin (f)) 9 

— j k J J N(r) | g — 2irR sinOdO rdr . (B-2) 


A-.B 

2ttR 


sin(j) 


Using the relationship sin0 = 2sin(|-) cos (|) , Equation (B-2) 
can be rewritten as: 


CO TT 

O s = 2A-^B N(r) J sin(2kr sin (|-) ) cos (|-) d0 rdr . 


(B-3) 


Letting u = 2kr sin(|) , Equation (B-3) can be expressed as : 
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(B -4) 


2A x Bk 3 jNCr) f sin(u) j^rdr 


2A,Bk 2 j N(r)[ 1-cos (2kr) ] dr. 
x o 

Asstaming N(r) = e ~l r//a ! , the net cross sectional area is 


a = 2A, B k 
s i 


J" e _r ^ a dr - J e _r/,a cos(2rk)dr 


(B-5) 


Performing the above integration and simplifying the result, 
the net cross sectional area is 


a = 2A,Bk a 
s 1 


l+(2ka) z J 


A B 4a 3 k^ 
l+(2ka) 2 


(B-6) 


where A in Equation (4-11-a) is 2A^ and B is [Ap 2 /p 2 + AB^/B 2 q ] 


According to Equation (B-5) the first term in the brackets 
is independent of frequency and a function of the average cor- 
relation distance "a", while the second term is a function of 
frequency (k=2rrf/c) and "a". It is the second term in Equation 
(B-5) that determines the attenuation dependency on frequency 
through parameter k. 
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FIGURE 2. The scattering configuration. 
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FIGURE 3. Plot of f-^(ka) and f 2 (ka) versus ka . 
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Z = per unit impedance 
Y = per unit admittance 
G = per unit conductance 


pAx is equivalent to per unit 
inductance 

Ax/B^ is equivalent to per unit 
capacitance 


FIGURE 4. A one-dimensional electrical model 

of a solid material. 
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FREQUENCY (JtHz) 

FIGURE 5. Attenuation versus frequency. 
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(b) Two realizations of mass density. 
FIGURE A-l. Model of a solid medium. 
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FIGURE A- 2 . Typical realization with the 
average and fluctuation in mass density 
versus distance. 
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